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Abstract 

In this note, we look at some hypoelliptic operators arising from nilpotent rank 
2 Lie algebras. In particular, we concentrate on the diffusion generated by three 
Brownian motions and their three Levy areas, which is the simplest extension of the 
Laplacian on the Heisenberg group H. In order to study contraction properties of the 
heat kernel, we show that, as in the case of the Heisenberg group, the restriction of 
the sub-Laplace operator acting on radial functions (which arc defined in some precise 
way in the core of the paper) satisfies a non negative Ricci curvature condition (more 
precisely a CD(0,oo) inequality), whereas the operator itself does not satisfy any 
CD(r, oo) inequality. From this we may deduce some useful, sharp gradient bounds 
for the associated heat kernel. 
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1 Introduction 

In the study of the long (or small) time behavior ( e.g. gradient estimates, ergodicity etc.) 
of simple linear parabolic evolution equations, one often uses lower bounds on the Ricci 
curvature associated to the generator of the heat kernel, see for example [T|llOt[TT] and the 
references therein. But this method fails in general in hypoelliptic evolution equations, 
since the Ricci (T2-) curvature in even the simplest example of the Heisenberg group can 
not be bounded below as explained e.g. in [9j[2]. Nevertheless, in the Heisenberg group 
case, many properties of the elliptic case remain true, and we shall details later some of 
the most interesting ones. 

Let us recall first some basic facts. 



The elliptic case 

Let M be a complete Riemannian manifold of dimension n and let C := A + V/i, where A 
is the Laplace-Beltrami operator. For t > 0, denote by Pt the heat semigroup generated 
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by C (that is formally P% = exp(i£)). For smooth enough function f,g, one defines (see 

r(f,g) = |V/| 2 = ±{£fg - fCg - gCf), 

T 2 (/, f) = \ {CT(f, f) - 2T(f, £/)) = |VV/| 2 + (Ric - VVh)(Vf, Vf). 

We have the following well-known proposition, see Proposition 3.3 in [I]. 
Proposition A. For every real the following are equivalent 

(i) . CD(p,oo) holds. That is T 2 (f,f) > pT(fJ). 

(ii) . For t > 0, T(P t f,P t f) < e- 2 ^P t (r(/,/)). 

(iii) . For t > 0, T(P t f,P t f) l 2 < e-t*P t (T(f,f)*). 

Moreover, in [7], Engoulatov obtained the following gradient estimates for the heat 
kernels in Riemannian manifolds. 

Theorem B. Let M be a complete Riemannian of dimension n with Ricci curvature 
bounded from below, Ric(M) > — p, p > 0. 

(i) . Suppose a non-collapsing condition is satisfies on M, namely, there exist to > 0, 

and z^o > 0, such that for any x £ M, the volume of the geodesic ball of radius to 
centered at x is not too small, Vol(B x (to)) > v$. Then there exist two constants 
C(p,n, vo, to) and C(to) > 0, such that 

| V log H(t, x,y)\ < C(p, n, u ,to) (^^- + ^ 

uniformly on (0, C(to)] x M x M, where d(x, y) is the Riemannian distance between 
x and y. 

(ii) . Suppose that M has a diameter bounded by D, Then there exists a constant C(p, n) 

such that 

\V log H(t,x,y)\ <C{p,n) (j + + ^ , 
uniformly on (0, oo) x M x M. 

The three-dimensional model groups 

In recent year, some focus has been set on some degenerate (hypoelliptic) situations, where 
the methods used for the elliptic case do not apply. Among the simplest examples of such 
situation are the three-dimensional groups (5 with Lie algebra q, where there is a basis 
{X, Y, Z} of such that 

[X, Y) = Z, [Z, Y] = aY, [Y, Z] = aX, (1.1) 
where a£l. The analysis reduces mainly to the thre cases a = 0, a = 1, a = —1. 

Example 1.1 (Heisenberg group, a = 0). The Heisenberg group can be seen the Euclidean 
space M 3 with a group structure o, which is defined, for x = (x,y,z),y = (x',y',z') £ M 3 , 
by 

xoy= ( x + x',y + y',z + z' + -(xy'-x'y) 
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The left invariant vector fields which are given by 
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x(f) = to - — 7 = ^ - |a. ) /, 

Y(f) = to J """"7" = ^ + f 8, ) /, 

£— >-0 



The right invariant ones are: 

_ /((e.0,0) of) -/(f) _ , + . 

£^0 £ V 2 

m = to —-"'7' = ^ - f a.) i, 

The Lie algebra structure is described by the identities [X, Y] = Z, [X, Z] = [Y, Z] = 0. 
In fact, all group structures satisfying U.l\) with a = can be transformed to the case 
(R 3 , o) by the exponential maps, the vectors fields {X, Y, Z} corresponding to the left ones, 
see Lemma 4-1 in JEj, see also /5]/- The natural sublaplacian operator for this model is 
C = X 2 + Y 2 . In this case, symmetries play an essential role : they are described by 
the Lie algebra of the vector fields that commute with C A basis of this Lie algebra is 
(X, Y, Z) and 9 = xd y — yd x . The last one reflects the rotational invariance of C, see JE/. 
For this sublaplacian C, we have 

T(f,f) = (Xf) 2 + (Yf) 2 , 

and 

r 2 (/, /) = (X 2 /) 2 + (Y 2 f) 2 + l(XYf + YXff + X -{Zf) 2 + 2(XZfYf - YZfXf). 

The appearance of the mixed term XZfYf—YZfXf prevents the existence of any constant 
pGi such that T2 > pT. Therefore the methods used in the elliptic case to prove gradient 
bounds cannot be used here. Nevertheless, B. Driver and T. Melcher proved in |6j/, the 
existence of a finite positive constant C2 such that 

V/ G C°°(e, M), Vt > 0, T(P t f, P t f) < C 2 P t T(f, /), (1.2) 

where Pt denotes the associated heat semigroup generated by C, C°°(]HI, R) is the class 
of smooth function form H to R with all partial derivatives of polynomial growth. More 
recently, H. Q. Li jlllf showed that there exists positive constant C\ such that 

v/ep°°(e), vt>o, T(p t f,p t f^ <CiP t (r(/,/)l). (1.3) 

(See also D. Bakry et al. JM/ for alternate proofs.) The gradient estimate \1.3\) is much 
stronger than 111. ffj) . and has many conseguence in terms of functional inequalities for 
the heat kernel Pt, including Poincare inequalities, Gross logarithmic Sobolev inequalities, 
Cheeger type inequalities, and Bobkov type inequalities, see section 6 in J||/. 

Let pt be the heat kernel of Pt at with respect to Lebesgue measures on R 3 . In 
H. Q. Li has also pointed out that for t > 0, g £ M, there exists a positive constant C such 
that 

IVlogptlfo) < (1.4) 

where d(g) denotes the Carnot-Cartheodory distance (see U.8\) ) between and g. This 
gradient estimate is sharp and plays an important role in the proof of hl.3\) . 
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In the case a = 1, the Lie algebra is the one of the SU(2) Lie group, and this case 
has been studied by F. Baudoin and M. Bonnefont in [4J. They show that a modified 
form of (|1.3j) and (|1-4|) hold. Other generalizations of Heisenberg group are the so-called 
Heisenberg type group. They have been studied by H. Q. Li in [12} I13j. where he shows 
that (|1.3|) and (|1.4|) hold in this setting. In this note, we shall focus on a group that we 
may call, the three Brownian motions model. It can be seen an another typical simpe 
example of hypoelliptic operator, but the structure is more complex than the Heisenberg 
(type) groups and the method of H.Q. Li fails to study the precise gradient bounds in this 
context. 

For this model, we shall first look at the symmetries, that is characterize all the vector 
fields which commute with the sublaplacian operator C, see Proposition ^. 11 The infinites- 
imal rotations are those vector fields which vanish at and a radial function is a function 
which vanishes on infinitesimal rotations. In this case, although the Ricci curvature is 
everywhere — oo, refer to [HIE], we shall prove that the T2 curvature is still positive along 
the radial directions, as it is the case for the Heisenberg group, see Proposition 13.11 As 
a consequence, the same form of gradient estimate (jl.4l) holds by combining the method 
developed by F. Baudoin and M. Bonnefont in [4] with the method in [12] . It is worth 
recalling that in [3j, D. Bakry et al. have obtained the Li-Yau type gradient estimates for 
the three dimensional model group by applying r2-techniques. In our setting, it is easy to 
see that this type of gradient estimate also holds. 

The three Brownian motions model 

The three Brownian motions model 913,2, see section 4 in [8], can be described as the 
Euclidean space R 6 with a the following group structure o, which is defined by for x = 

(x 1 ,x 2 ,x 3 ,y 1 ,y 2 ,y 3 ),y = (x l ,x 2 ,x 3 ,y[,y' 2 ,y 3 ) G M 6 , 

(xi,x 2 ,x 3 ,yi, 2/2,2/3) (4, 4,^3, 2/1,2/2, 2/3) = { x i +x'i,x 2 + x' 2 ,x 3 + x' 3 ,y 1 + y[ + -{x 2 x 3 - x 3 x' 2 ), 

V2 + y' 2 + ^{xzx' x - xix' 3 ),y 3 + y 3 + ^{xix' 2 - x 2 x' x )). 

For simplification, we make the convention that the index i = j mod 3, and here we choose 
j = 1, 2, 3. In what follows, denote ^3,2 = (R 6 , °) be the three Brownian motions model. 
The three left invariant vector fields which are given, for 1 < i < 3, by 

f = nm /(fo( £l , £2 , £3 ,0,0,0))-/(f) = / _ x_^ ^ + Xj ^ bi x 
= nm f(xo(0,0Ae 1 ,£ 2 ,£ 3 ))-f(x) = B 

s^O £ 

where £% = £ and £j = for j ^ i. Here we use the notation d{ = d Vi . 

The right invariant vector fields which are gives JQ, for 1 < i < 3, £% = £ and £j = for 

j + i, 

± = nm /(( £l , £2 , £ 3,0,0,0)of)-/(f) = , + _ X _^ B , 

£-!>0 £ \ 2 2 / 

There are no Y^s since in this setting the left and right multiplications coincide. The Lie 
algebra structure is described by the formulae, for 1 < i,j < 3, 

[Xi,X i+1 ]=Y i+2 , [X i ,Y j } = 0. (1.5) 
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Similarly for all group structure satisfying (jl.5p can be transformed to the case (M 6 , o) via 
the exponential maps, the vectors fields are corresponding to the left ones. 

In what follows, we are interested in the natural sublaplacian for this model, which is 
defined by 

£ = t X l 

i=l 

The reason why we call it the three Brownian motions model is that is the infinitesi- 
mal generator of the Markov process ({-Bi}i<«<3' {\ J * BidBi + \ — £?j + idBj}i<j<3), where 
{-B,;}i<i<3 are three real standard independent Brownian motions. 

For all t > 0, P t := e tc denotes the associated heat semigroup generated by £,, pt the 
heat kernel of Pt at with respect to Lebesgue measures on M 6 . For this operator C, we 
have 

3 
i=l 

and 

3 3 

r 2 (/,/) = Yl ( x * x jf) 2 -2Y,Xif(X i+1 Y t+2 f - Y i+1 X i+2 f). 

i,j=l i=l 

Here again the mixed term Yli=i — ^+1-^+2/) prevents the existence of 

any constant p such that the curvature dimensional condition CD(p, 00) holds. Neverthe- 
less, we have the following Driver- Melcher inequality, see |15| . 

T(P t f,P t f)<CP t (T(fJ)), 

for some positive constant C. The constant C here can be expressed explicitly following 
the method in [2]. Also the optimal reverse local Poincare inequality holds, see Remark 
3.3 in [2J. That is, for any t > and any / E C™{m 3 , 2 ), 

tT(P t f,P t f)<^(P t (f)-(P t f) 2 ). 

For the H. Q. Li inequality (|1 .3[) . the methods deeply rely on the precise estimates on 
the heat kernel pt and its differentials (see [111 [2]). Up to the author's knowledge, these 
precise estimates are not known in the three Brownian motions model, neither the H. Q. 
Li inequality. Nevertheless, we shall prove that one of the key gradient estimates p. 41) 
holds, which would be a first step for the proof of the H. Q. Li inequality in this context, 
see Proposition 14.21 

The dilation operator in this model is defined by V := | Yl^ X A + Yli=i vSii an d ^ 
satisfies 

[C,V]=C. (1.6) 
For t > 0, let T t = e tv be the semigroup generated by T>, that is 

T t f (xi,x 2 , x 3> yi, 3/2. 3/3) - / 

From the commutaton relation (jl.6|) . one deduces, for t, s > 0, 

PtT s = T s P e s t . 
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Since is a fixed point of the dilation group Tt, it follows 

P t (/)(0) = P 1 (T logf /)(0). (1.7) 

So it is enough to describe the heat kernel at any tme and any point to know the operator 
A(/)(0). 

The natural distance, induced by the sublaplacian operator C, is the Carnot-Caratheodory 
distance d. As usual, it can be defined from the operator C only by 

d(gi,92)-= sup f(gi) - f(g 2 ). (1.8) 
{/=r(/)<i} 

For this distance, we have the invariant and scaling properties, see [HI I17j. 

d{gx,gz) = d^ 1 05-1,0) := d(g 2 l o #1), and d(^x^ 2 y) = jd(x,y), 
for all gi,g 2 G 9T 3i2 , 7 G R + and x = (x 1 ,x 2 ,x 3 ),y = (1/1,3/2,3/3) G 

2 Rotation vectors and Radial functions 

In this section, we shall characterize all the vector fields which commute with C Obviously 
the right invariant vector fields {Xi, ii}j=i 2 3 commute with C since they commute with 
{Xj}i=i 2 3. Like the rotation vector field = xd,, — yd x in the Heisenberg group, which 
commutes with C, there are three rotation vector fields in this case 

9i = x i+ id i+2 - x i+2 d i+ i + yi+id i+ 2 - y i+2 d i+ i, i = 1, 2, 3. 

It is easy to see that {#i}i=i, 2 ,3 commute with £ and we have [9i, 0i+i] = @i+2, for 1 < i < 3. 
We first have the 

Proposition 2.1. T/ie vector fields which commute with C are the linear combination 
of the following nine vector fields: the three right invariant vectors, the three rotations 
{#i}i<i<3, and d\, d 2 , d 3 , that is 

T = {X : X G span{Xi,Yi, 1 < i < 3}, [£, X] = 0} = Linear{X, 0;, Yj, 1 < i < 3}. 

(2.1) 

Here "span" means the we consider linear combinations of the vector fields with smoth 
functions as coefficients, , while "Linear" means the coefficients are constants. 

Proof. We only need to show that the left hand side space in (I2.ip is contained in the right 
hand side one. To this end, for any vector field X = X^i=i a iXi + hY^ for some smooth 
function a,, hi, satisfies [£, X] = 0. For 1 < i < 3, denote Z$ = XiXi + \ +Xi + \Xi, it yields 
X,- L X l+l = Z ' +2 + y '+ 2 and XiX i+2 = Z ' +1 ~ y ' +1 . Notice that 

3 

[C,X] = ^2 (y£- a iXi + (Cbi + — Xj +2 ai + i)yj + (Xi + iai +2 + Xj +2 aj + i)Zj 

+ 2X iai Xf + 2X i b i X i Y i + 2(Xib i+1 - a l+2 )XiY i+1 + 2(XA+2 + a m )A^+ 2 ) , 
thus we have, for 1 < i < 3, 

A>; = = 0, (2.2) 
Xi+\bi + 2 = —Xi +2 bi + i = ai, (12.21) 
XiOi+i = —Xi+iai. (|2.2I ; ) 

Let us first prove the following two claims. 
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Claim I: For 1 < i < 3, aj is independent on {yi, 1 < i < 3} and linear in {xi, 1 < i < 3}. 
To proof the desired result, for 1 < i < 3, we have the following commutative 
property: pQ,Y i+1 ] = 0, together with (|2.2|) . it yields XjYj +1 6j = 0. Since Yj + i = 
[Xj + 2, X], we can get XfX i+2 bi = 0, thus Xfa i+ i = by the relation X i+2 bi = a i+ i. 
Similarly we have Xfdi +2 = 0. In fine, together with XjOj = 0, 



Since 



we have 



Xf aj = 0, i,j = 1,2,3. 



[Xi,Y 3 ]b 2 = 2X!X 2 a 3 = -X 2 X 3 ai, 
[X 3 , Yi]6 2 = 2X 3 X ia2 = -X 2 X 3 ai, 
[Xi,y 2 ]6 3 = 2X!X 2 a 3 = -X 3 X ia2 , 

XiX 2 a 3 = X 2 X 3 ai = X 3 Xia 2 = 0. 



Together with the fact Xidj = —Xjdi by A2.21 7 ) . we have XiXjdt = for i,j,k all 
different. Thus we can conclude 

XiXjdk = for 1 < k < 3. (2.3) 

Note that Ylai = X 2 X 3 ai — X 3 X 2 ai, Y 2 ai = Xfa^ and Y 3 ai = —X\d 2 , thanks to 
(|2.3p . we get Y^ai = 0,i = 1,2,3. That is a\ is independent on {yi,i = 1,2,3}. 
Similarly a 2 ,a 3 is independent of {yi,i = 1,2,3}. Then from the definition Xi, we 
have Xidj = didj for 1 < i, j < 3. With (|2.3p . we can conclude that {cij,i = 1,2,3} 
is linear in {xi,i = 1, 2, 3}. 

Note that we can also write in the form X = Yli=i a i^i + c $i-, where Ci = b{ + 
|(ai_|_ 2 Xj + i — aj + iXj +2 ), 1 < i < 3. Then we can conclude 

Claim II: {cj,z = 1,2,3} is linear in {xi,yi, 1 < i < 3}. By Claim I, a, is independent 
of yj, together with the fact Xidi = 0, we have the equation (|2.2I ) is equivalent to 

-Oi = X + iCj +2 + -Xi + \di + \ai 

2 2 ! (2-4) 

= -x l+2Cl+1 + -x i+2 a J+2ai . 

And Xibi = is equivalent to 

XjC, = ^(xi + i5ia i+2 - Xi +2 dia i+ i). (2.5) 

Using ()2.2[) - ()2.5|) . the relations [Xj, Xj+i] = Yi +2 and Claim I, through computation, 
we have 

YiCj = didj, for i,j = 1,2,3. (2-6) 

Since di is linear in Xi, we can conclude that Cj has no second order terms in 
{xi, yi, 1 < i < 3}. By the definition of Xi and (|2.5p and (|2.6p . we have 

diCi = Xid H —Y i+2 Ci — Yi+iCi 

= 2(^+1^^+2 - x i+2 did i+ i) - ^Y-did i+2 + ^Y-did i+ i 
= 0. 
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similarly, 



By (Ha]) and UM, 

OiCi+i — AjCj + i H — Yi + 2Ci + i — Yi + iCi + i 

1 1 1 ( 2 - 7 ) 

— 2 ai+2 ~ 2 Xi ® iCLi+2 ~ 2 Xi+l ® i+iai+21 

diC i+2 = --<H+\ + 2 X i^°j+i + 2 Xi + 2di + 2ai + 1 - ( 2-8 ) 

By Claim I, ()2.6|) - (|2.8j) and diOi = 0, we can conclude that for 1 < i,j < 3, c^cj is 
constant. Thus we complete to proof Claim II. 

By the above two claims and <9jOj =0, we can assume 

di = Ai y i + \Xi + \ + Ai y i + 2Xi + 2 + Bi, 
where A% j = —Aji, Bi are constants, then we have, by (|2.6p - (]2.8p . 

Ci = -{B i+ iXi + 2 — B i+ 2X i+ i) + A i:i+ iy i+ i — A i+2 ,iyi+2 + Di, 

where Di are constants. 

If we choose Bi = 1 (or respectively Di = 1, ^^+2 = 1) and the other constants 0, we 
get X = Xi (or respectively Yi, 9i). Thus we complete the proof. 

□ 

In the Heisenberg group, the radial functions / can be characterized by 9f = 0. Here 
in our setting, as an extension of such characterization, we can give a definition of radial 
functions. 

Definition 2.2. A smooth enough function f is called radial if and only if for 1 < i < 3, 
6if = 0. 

(Notice that here the vector fields 9i are the commuting vector fields which vanish in 

0.) 

Remarks 2.3. Note that the heat kernel (pt)t>o is radial. The reason is that for any 
function f, 1 < i < 3, 6if(0) = and {0i}i<i<3 commute with C, whence they commute 
with the semigroup Pt = e . Hence, for any function f , one has PtOif = 0, which, taking 
the adjoint of Bi under the Lebesgue measure, which is —di, shows that for the density pt 
of the heat kernel at 0, one has 9iPt = 0. This explains why any information about the 
radial functions in turns give information on the heat kernel itself. 

Remarks 2.4. For any radial function f, there exist some function g such that f(x, y) = 
g(ri,r 2 ,z), where n = X^Li^ 2 ' r 2 = YA=iVh z = Y?i=\ x iVi- Indeed, by the definition, 
for 1 < i < 3, 9if = 0, then we have f = f(Ux, Uy), where U is arbitrary linear orthogonal 
transformation on R 3 ; which satisfying U*U = UU* = 1. Hence f = f{\x\, |y|, (x,y)), for 
some function f. Here is another way, by the transformation 9i, we will directly get 

z 



f(x,y) =/(Vn,0,0,— ,0,Vr2-z 2 /n). 

vn 



s 



3 T2 curvature 

Recall that we can't find a constant p £ R such that T2 > pY because of the appearance of 
the items XifXjY^f . In other words, the Ricci curvature is everywhere —00. Nevertheless 
we shall prove I^ curvature is positive on the radial functions. 

Proposition 3.1. For any smooth radial function f, we have 

r 2 (/,/)>o. 

Here we will give two different proofs. The first one is that we shall use directly 
the three equations asserting that a function is radial. Then, applying the vector fields 
{Xj}i<j<3 on these equations, we get nine equations in hand. It follows that we can get 
the exact expressions of {^i^/}i<i,j<3 in terms of XiXjf and also first order terms. (In 
fact, we adapt the mathematical software MAPLE to do it). Then we substitute them 
into the formal expression of I^, and we find that Y2 can also be expressed in a functional 
non negative quadratic form. 

The second way is that by the Remark 12.41 we have an expression of the sublaplacian 
operator acting on radial functions directly through a good parametrization, say r\,r-i,z. 
Through computation, we can obtain the exact expression of T2 curvature and find again 
that T2 can be expressed in a functional non negative quadratic form, thus we are done. 

The first proof. A radial function / satisfies 9if = 0, which is equivalent to say that 

v f Y , x i+l + x i+2 v , , XiXi+l ~tyi+2 v XiX i+2 + 2y i+1 
x i+1 X i+ 2f-x i+ 2X i+1 f ^_ — ?-Yif-\ - Y i+2 f = 0. 

Differentiating the above equations at the directions {Xj}i<j<3, with the commutative 
relations = Y i+2 , we get the nine differential equations, for 1 < i < 3, 

^+1^+2^1/ — Xi + 2Xi + \Xif XiYif 

. 2^Ei+l-2yi+2 XjX i+ 2 + 2y i+ i 
H 2 XiY i+1 f H A1Y3/ = U, 

Xi+2f + Xi + iXi +2 Xi + if - x i+2 X i+1 f X i+ {Yif 

XiX i+1 - 2y i+2 XiX i+2 + 2y i+1 
H 5 X i+1 Y i+1 J H X i+ iY i+ 2j = U, 

— X i+ \f + Xi + \X i+2 f — x i+ 2X i+ iX i+ 2f X i+2 Yif 

XjXj+i - 2y i+2 XjXj+2 + 2y i+ i 
H 5 X i+ 2Yi + \f H Aj +2 Yj+2/ = U. 

For simplificity, we will use the following notations, for 1 < i < 3, 

a i+ i := x i+ iXif - XiXi+if, f3 i+ \ := yi+iXif - yiX i+i f, 

li ■= Xi,y i+ i - Xi+iyi, rji := xiyi + x i+ iy i+ i, 
\x\f := Xi + xf +1 , Ai := ~f 2 XiX i+ if + ^X 2 X i+ if + 71X3X^+1/, 
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and 

3 

\h\ 2 := N 2 > for h = (^'^2, h 3 ) € M 3 . 
From the above nine differential equations, we can get, for 1 < i < 3, 

XiY i+1 f = -^-^(^(xiX i+ i\x\ 2 + 2y i+2 \x\ 2 - 2x i+2 rn + ^ViVi+i) ■ (x i+2 Xi +l Xif - x i+ iX i+2 Xif) 

+ {x i+1 x i+2 \x\ 2 - 2yi\x\ 2 +1 + 2xi7] i+1 + Ay i+ iy i+2 ) ■ (x i+1 X 2 f - XiX i+1 Xif - X i+1 f) 
+ (x 2 +1 \x\ 2 + 4y 2 +1 ) ■ (xiX i+2 Xif - x i+2 X 2 f + X i+2 f)) , 

and 

XiY i+2 f = -^-^{^XiX i+2 \x\ 2 - 2y i+1 \x\ 2 +2 + 2x i+ i7] i+2 + % +2 y;) • {x i+2 X i+ iXif - x i+ iX i+2 Xif) 

+ (x i+1 x i+2 \x\ 2 + 2yi\x\ 2 +1 - 2xiT] i+ i + Ay i+1 y i+2 ) ■ (xiX i+2 Xif - x i+2 X 2 f + X i+2 f) 
+ (x 2 +2 \x\ 2 + 4y 2 +2 ) • (xi+iXff - Xi X i+l Xif - X i+1 fj) . 

Note that 

3 3 

r 2 (/, /) = (Wjf) 2 -2j2 X if(X i+1 Y i+2 f - Y i+1 X i+2 f). 

i,j=l i=l 

With the exact expressions of XjYjf in hand, through calculation, we have 

3 

F : = 2| 7 | 2 ■Y J Xif{X i+1 Y i+2 f - X l+2 Y l+1 f) 
i=i 

3 

= 2 Y (Xi+jf(2y i+j+1 ii +1 - x i+j+l yi\x\ 2 i+1 + XiX i+j+1 r]i +1 ) 

- X i+j+1 f(2y i+j 7 i+1 - x i+j yi\x\ 2 +l + XiXi +j r]i +1 )^XiX i+j+2 f - \x\ 2 ■ \a\ 2 - 4|/3| 2 

3 

= 2 ^ (2f3 i+j -f i+1 + a i+j jiX i+1 - ai +j ji +2 x i+2 ^XiX i+j+1 f - \x\ 2 ■ \a\ 2 - 4|/3| 2 . 
Rearrange the items, we have 

3 3 

F = 2 Y aiX i+j+1 (~f i+j X i+j X i+1 f - li+j+l X i+j+2 X i+l f) + 4^/?^ - \x\ 2 ■ \a\ 2 - 4|/3| 2 . 

i,j=l i=l 

Notice that 

3 3 3 

ItI 2 • ^2(XiXjf) 2 -^A 2 = ^2 (n +j X i+j X i+1 f - n +j+1 X i+j+2 X i+1 f) , 

i,j=l i=l «J=1 
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it follows that 

l7l 2 -r 2 = | 7 | 2 - j^iX^ff-F 

3 3 
= X^^ 2 ^ ~~ + ^ hi+3 X i+3 X i+tf ~ 7i+j+lXi+j +2 X i+1 f) 

i=l i,j=l 
3 

- 2 y~] aiXi + j + i(^i + jX i+ jX i+1 f - ji + j + iXi + j +2 X i+ if) + |x| 2 • |a| 2 

3 3 

= ^( 2 A - ^4j) 2 + ^ (7i+jX i+ jXj + i/ - ji + j + iXi + j +2 Xi + if - aix i+j+ i) 2 . 

i=l i,j=l 

Hence we complete the proof. 

□ 

The second proof . Denote r\ = Ylf=i x h r 2 = Y^i=iyh z = Ylf=i x iUi- Through calcu- 
lation, we have 

Cr\ = 6, Cr 2 = r±, Cz = 0, 

and 

r(n,ri) = 4n, r(r 2 ,r 2 ) = nr 2 - z 2 , T(z,z)=r 2 , 

r(n, r 2 ) = o, r(n, z) = 2z, r(r 2 , *) = o. 

For any radial functions f,g depend only on r±,r 2 ,z, we have 

Cf(n, r 2 , z) = drJCn + d r J£r 2 + dJCz + d 2 riri fT( ri , n) + d 2 2T . 2 fr(r 2 , r 2 ) + 8% z fT(z, z) 
+ 2d 2 riZ fT(n,z) + 2d 2 2 JT(r 2 , z) + 2d 2 ir JT( ri , r 2 ) 
= A n d 2 . ir J + (nr 2 - z 2 )d 2 2r J + r 2 8 zz f + 4zd 2 lZ f + 6d r J + n d r J 
■■=Cf, 

where £ has the following expression 

tf = 4n/n + (nr 2 - z 2 )/ 22 + r 2 f zz + 4z/i a + 6/i + n/ 2 . 
Hence for any radial function / = f(r±,r 2 ,z),g = g(ri,r 2 , z), 

T(f,g) := l -(C{fg)-fCg-gCf) 
= \{C{fg)- fCg-gCf) 



and also 



t(f,g), 



T 2 (/,/) :=^(CT(f,f)-2T(f,Cf)) 

= i(£f(/,/)-2f (/,£/)) 
= f 2 (/,/). 

11 



Through direct calculation, we have 

r(f,g)=f(f,g) 

= 4ri/i0i + (nr 2 - z 2 )f 2 g 2 + r 2 f z g z + 2zfig z + 2z/ z #i 

and 

r 2 (/,/) = f 2 (/,/) 

= Wrjfa + 16ri/i/u + 8n(nr 2 " ^)/i2 + 8 ( r i r 2 - ^)/ 2 /i2 + 8(nr 2 + z 2 )/^ 
+ 32-nzfufu + n(nr 2 - z 2 )/ 2 / 22 + (nr 2 - z 2 ) 2 / 2 2 + 2(nr 2 - z 2 )/ 2 / 22 

2 

+ 8z(nr 2 - z 2 )/ 12 / 2 . + (2r 2 + ^)/f + 2r 2 (nr 2 - z 2 )/2* + r 2 /J, + 4r 2 /i/„ 

+ 8r 2 zf lz f zz + 16zhf lz + 8z 2 f u fzz + 12/f + ^n/J - 4(nr 2 - z 2 )/i/ 22 
- 4n/i/ a - (rir 2 - z 2 )f 2 f zz - 2zf 2 f z . 
By careful study, we can express the above into a functional quadratic form. 

r 2 (/, /) = ((nra - z 2 )/ 22 + ^/ 2 - 2/i) 2 + 8n(nr 2 - z 2 ) (f 12 + A + ^A) 2 



2 ^ 2 N _c r l j. j. \ 2 / r l j. r,j- r l r 2 ~~ z 



+ - {(nr 2 - z*)f 2z + -ff z -zj 2 ) + [j h ~ 2/i - f, 



2 \ 2 



+ 4 ( A + ^/ 2Z + 2n/n + 2zf u I + 2(rir 2 - z 2 ) ( -/„ + 2/ la 



y 2 \ 2 

„— /« + 2ri/n + 2z/i J ...... ... . 

2ri / \ri 

Hence the desire result follows. □ 

4 Gradient bounds for the heat kernels 

As done in [3], we have the following Li-Yau type inequality holds. 

Proposition 4.1. There exist positive constants Ci,C 2 ,C^ such that for any positive 
function f , if u = log Ptf , we have 

3 r 
d t u > dT{u) + C 2 t J2 \YM 2 ~ -j-- 
i=i 

Proof. Here we briefly proof it for the readers' convenience. Notice that for all A > 0, 

3,3 3 3 



r 2 (/, /) = Y^(xff) 2 + 1 ^T(Yif) 2 + 2 (A, l+ i(/)) 2 + 2 £ (XifXi+iY^f - x i+2 /x,y i+1 /) 

i=l i=l i=l 

>\{cf) 2 + \jz^f?-^Vnf)- 



3 



\ i=i 



i=l i=l 
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where Di i+i = + and the last two inequalities follow from the Cauchy- 

Schwartz inequality. 

Set f s = P t - S f,u s = log/.s, following [3], let 

*i00 = p s (/ s rK,0) ; * 2 («) = ^ ^/»£(*i«,) 2 J , 

we have 

*!(«) = 2P s (/ s r 2 ( Us ,n s )), $ 2 ( s ) = 2P S ^/ s ^r(y^)^ . 

Combining (|4.ip and 

\2 \ o_,/*„. ..2 /v. ^/s 



(£u s ) > 27>Cn s - 7 , £u s = — T(u s ), 

Is 



we have 



Denote a, 6 are positive functions defined on [0, t), with b is decreasing, we have 



A 3 ; WM ' v ' v ' z 3 J 3 
By choosing 

, b 3b" 3b' 

a = - b > x = ~v> 7= 4^ + T' 

and then choose 6(s) = (i— s) a , for some a > 2, integrating the above differential inequality 
from to t, the desired result follows. □ 

As a consequence, we have the following Harnack inequality: There exist positive 
constants A\, A%, for i 2 > t\ > 0, and <?i,<7 2 G ^3,2, 

^^)<f^ Al /^. (4 . 2) 



Pt 2 {92) \t\ 

Here is an analogue result of Theorem B in the three Brownian motions model. 
Proposition 4.2. There exists a constant C > such that for t > 0, g = (x,y) G 9^3,2? 

vT(io gft )0,) < 

where pt(g) denotes the density of P t at and denotes the Carnot-Caratheodory 
distance between and g. 

Proof. Following [3], for < s < t, let <3?(s) = P s (pt-s log pt-s) j we have 

$'(«) = p s (pi_ s r(iog P i_ s )), $"( s ) = 2P s (p t _ s r 2 (iog Pt _ s )). 
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By Proposition 13.11 3>" is positive, whence $' is non-desceasing, thus 

/ &(s)ds > V(0). 
Jo 2 

That is 

2 

p t T(\ogp t ) < -{P t /2(Pt/2^ogp t / 2 ) -ptlogpt). 

The right hand side can be bounded by applying the above Harnack inequality (|4.2p and 
the basic fact Pt/2(d) — Pt/2^)-, f° r an 9 £ 9^3,2- We have 



v / r(iogp,)(5)<c(^ + i= 



In particular, 



VT(logp 1 )(g)<C(d(g) + l). 

If <^(<?) > 1, it is trivial to get the desired result. 

Note that for x = (x±, ^3), y = (j/i, 2/2,2/3), see P. 125, Theorem 1 in |8], 

Pi(x,y) = -(2vr)~^ ^ exp(-iy • a)-y- fsinhiy-J 

Xt4.^X^ / led I \ 1 t~t 

__ coth T j| JJ_ da k 



(*) /" ^ I ' X ' 

(27r) 2 / exp (— • a) — I smh ■ 
• expi |-|x| 2 + (|x| 2 - • d) 2 ) ^1 - ^-cothiyj I da k , 



2 r 

3 



where a = (ai,a 2 ,a 3 ), \a\ 2 = J2l=l a l> ® = y^|( a i' a 2, 03) and A = f 



k=l 

Ol — CL'Z 
-OL\ O3 

Q!2 — «3 

follows from x^4 2 x* = — |a| 2 (|x| 2 — [x ■ a) 2 ). Notice 

15 F f I \ ^ { Ictl 

djpi(x,y) = (2-7r)~"2" / exp(—iy ■ a)\a\ sinh — • — a%- + [xa — x ■ dd,)(l coth 

Jm? V 2 / V 2 

• expi |-|x| 2 + (|x| 2 - (x • a) 2 ) ^1 - ^y-cothy^ j da* 



and 



%>i(:c, y) = (2vr) 2% ^ ^ exp • a)^ ^sinh ^\ ■ at 

■ exp^ |-|a;| 2 + (|x| 2 - (x ■ a) 2 ) ^1 - coth-y-^ | da fc . 

Let 

Wl= [ \a\ (sinh^- ) TT da fe , W 2 = [ (sinh^-J • \a\ 2 coth TT cfa fc , 

* \ 2 / fc =1 JR3 V 2 / 2 ^ 
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obviously we have W\,W 2 are bounded. For g = (x,y) £ 9X3,2, satisfying d{g) < 1, by the 
basic fact that \x\ < d(g) < 1 (see [HJ[T7], i n fact we can easily proof it on the nilpotent 
groups.), we have 

vT(p)(s) < C x \x\(Wx + W 2 ) < C 2 \x\ < Cd{g). 
The desired result follows by the time scaling property (| 1 . 7j) . 

□ 

Notice that for any radial function /, Ptf is also radial since all 6i commute with 
P%. Thanks to Proposition 13. R we have H. Q. Li inequality, LSI inequality, isoperimetric 
inequalities etc. hold for the semigroup restricted on the radial functions, see pOU]. We 
state them in the following proposition. 

Proposition 4.3. For any compactly supported smooth, radial function f , for any t > 0, 
9 G 9*3,2, 

(%) H. Q. Li inequality. r(P t f,P t f)*(g) < P t (T(f,f)*)(g). 

(ii) LSI inequality. P t (f log f)(g) - P t {f) log P t (f)(g) < tP t (^) (g). 

(Hi) Isoperimetric inequality. P t (\f - P t (f)(g)\)(g) < Ay/iP t (T{f)^){g). 

Discussion: Here we have shown that H. Q. Li inequality holds for the radial func- 
tions. For the general functions, it is still open. Following the viewpoint of |llj . also [2], 
one key point to proof H. Q. Li inequality is the precise lower and upper bounds for the 
associated heat kernel. But in our setting, this estimates is unknown, at least the methods 
in |llj-|13j are not applicable. This precise estimates are also essential to proof the cheeger 
type inequality, see Lemma 5.1 in [2] and proof the constant coefficient is bounded in the 
complex quasi-communication method, see Proposition 5.5 in [2]. 
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